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1 Introduction 

Since this is the first contribution to this meeting which is mainly involved in cosmology, I 
would like to put it into its cosmological perspective: 

In the standard model of cosmology we assume that the universe is homogeneous and 
isotropic on large enough scales. This assumption, originally probably made just by reasons 
of simplicity, leads to a so called Friedmann Lemaitre universe which explains, e.g., the well 
established uniform Hubble expansion |]J. Within this model, the universe started out from an 
extremely dense, hot initial state and subsequently cooled by adiabatic expansion undergoing a 
series of phase transitions. At a temperature of about 0.1MeV« 10 9 K, deuterons become stable 
and virtually all the neutrons present are bound into 4 He. This leads to the well established 
abundances of the light elements. As the universe cools further, below about 3000°K there are 
no longer enough ionizing photons around to keep the Hydrogen Helium plasma ionized. The 
matter in the universe recombines to neutral H/He and the universe becomes transparent to 
the cosmic radiation field. Radiation can then propagate freely, influenced only by the cosmic 
expansion and redshifted to the 2.7°K microwave background which we observe today. At re- 
combination, the age of the universe was t ~ 2 x 10 5 years, which corresponds to a redshift of 
z R « 1100. 

We now want to discuss the possibility, that the cosmic plasma might have reionized again 
at some lower redshift Zj < zr. The reason one might want to investigate this idea is twofold: 

• Observationally: The Gunn Peterson test, i.e. the absence of a Lyman-a trough in all 
the observed quasar spectra, shows that the intergalactic medium is ionized for z < 4.5 




• Theoretically: The anisotropies in the cosmic microwave background (CMB) have 
turned out to represent one of the most stringent 'bottle necks' which a scenario of 
structure formation has to pass in order to be acceptable. As a possibility to relax this 
constraint, it has been proposed that early reionization can damp CMB fluctuations on 
small scales (6 < 6°) due to photon diffusion in the ionized plasma We shall illustrate 
this idea in this paper. 

In the next section, we want to explain the influences of reionization on the cosmic microwave 
background: It leads to damping of fluctuations due to photon diffusion on one hand and slightly 
distorts the CMB spectrum on the other hand. In Section 3, we discuss the processes in the 
plasma which have to be taken into account to describe reionization and formulate the system 
of differential equations. Finally we draw conclusions and outline future progress. 

In this contribution we use the following notation: 

• Greek indices run from to 3 and latin ones from 1 to 3 



• We assume throughout a spatially flat universe, i.e. the density parameter, Qtot = 
Ptot/Pc = 1 (a possible curvature would only change angles, but no other conclusions 
since we are mainly interested in the regime z » 1). 

• We choose t as conformal time coordinate and work with the metric signature (—,+,+,+), 
so that ds 2 = a(t) 2 (-dt 2 + dx 2 ). 

• Boldface characters denote 3 dimensional vectors. 

• We parametrize Hubble's constant by Hq = 100km/(s Mpc)/i. The density parameter of 
the baryons, is denoted by Qb = Pb/Pc = Pb/(7-7^5o x W~ 29 g / cm 3 ) . Observations limit 
0.5 < h < 0.8 and (including nucleosythesis calculations) 0.01 < £Ib < 0.1. 



2 Reionization and the cosmic microwave background 
2.1 Damping of fluctuations by photon diffusion 

As we shall see in this paragraph, reionization can lead to a substantial damping of fluctua- 
tions in the CMB on small angular scales. This is important for some scenarios of structure 
formation to overcome present limits posed by small and medium angular scale experiments 
(see contribution of P. Richards in this proceedings). As an example, we mention the recently 
investigated scenario with cold dark matter (CDM) and texture seeds || |(| |?]) H H- There, 
an analysis of CMB anisotropies shows that early reionization is a crucial ingredient for this 
scenario [^0[. Without damping, the small scale anisotropies would dominate and exceed ob- 
served limits. On the other hand, studies of the texture scenario show, that textures lead to 
early formation of objects (7j. At a redshift z ~ 50 about 1% of the baryons in the universe 



have collapsed and formed objects with mass M n i < 1O 5 M0 [11]. If one assumes that, due to 
the formation of these objects, radiation energy of about lOOkeV per nucleon is emitted, this 
would yield a total energy density of this ionizing radiation p. L = qps ~ 10 -6 /9b ~ 10 3 eV x ng, 
which is by far enough to reionize all the hydrogen in the universe. 

By early reionization in this context we mean that it has to happen early enough so that 
Compton scattering of electrons is still effective. At late times electrons are too sparse to 
scatter photons effectively and the electron proton plasma is again decoupled from radiation. 
This decoupling time is determined by the optical depth due to Compton scattering being equal 
to unity: 

f^dec 

T = / n e <JTadt = 1 , 
Jt 

which leads to a decoupling redshift of 

where aj- denotes the Thomson cross section. We further assume (for simplicity) £l tot = 1, so 
that a substantial amount (about 90-95%) of non baryonic, dark matter has to be present. We 
parametrize the metric of this spatially flat Friedmann universe by using conformal time, 

ds 2 = a 2 (-dt 2 + dx 2 ) . (2) 

The physical time differential is thus given by adt. 

We now want to investigate how CMB anisotropies on scales smaller than t^ ec , i.e., 6 < 
1 / \J Zdec + 1 ~ 6° can be damped if reionization happens substantially before Zd ec - 
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^/-r>v^- = c[/] , (3) 



Boltzmann's equation for Compton scattering is given by 

df_ 

dp 1 

where / denotes the distribution function of the photons and C[f] is the collision integral. 
In a perturbed Friedmann universe we separate / into an isotropic background contribution 
and a small perturbation / = / + Sf. One can then find a gauge-invariant (i.e. invariant 
under linearized coordinate transformations) variable which reduces to 5f for perturbations 
which are much smaller than the size of the horizon. For the energy integrated "brightness 
perturbation" 

47T 



M 

P 



roc 

/ Tp^dp 
Jo 

the perturbation of the Boltzmann equation (||) then becomes (for scalar perturbations) 



M + e*diM = 4e i a i ($ - *) + aa T n e [D^ -M- AeHdiV + -&.-AT J "] . (4) 

y 2 

Here e is a unit vector which denotes the direction of the photon momentum, V is a potential 
for the baryon velocity, <£, are the so called Bardeen potentials which parametrize scalar 
perturbations of the geometry and I is an arbitrary length scale introduced merely to keep the 
perturbation variables dimensionless. Furthermore 

-lx 

flW = (l/4vr) J M(e)dn and 

Mij = A f M{e)e ij dn] . 
The baryon equation of motion is 

IdjV + (a/a)WiV = di$ - aaTriePr (Mj + Ald { V) , (5) 

3pb 
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with Mj = — ejMdQ, 



An introduction to gauge-invariant cosmological perturbation theory and a thorough deriva- 
tion of eqs. (Q,||) is given in Q. 

In eq. (Q) the first contribution to the right hand side is the gravitational force. In the 
second contribution, the collision integral, the first two terms are the usual smoothing, then 
there is a Doppler term and the final term is due to the anisotropy of the Compton cross section 
(it can be neglected in the limit of very many collisions). In eq. (||) we see, that on the right 
hand side the drag force due to the radiation drag experienced by the electrons is added to the 
well known gravitational acceleration term. 

We now want to estimate the damping of M due to photon diffusion. Setting tx = 
(ao"T^ e ) _1 ) w e investigate the limit tx/t « 1. In deepest order we find 



I>W = (4/3)lAV = (4/3)Z)<, s ) , 
second equal sign is due to the baryon number conservation equation. In this limit therefore, 



where Dg is a perturbation variable for the baryon energy density perturbation, and the 
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the entropy per baryon is conserved, 5s/ s = Sub/tib, i-e., baryons and radiation behave like 
a single perfect fluid with energy density p r + ps and pressure p r = (l/3)p r . We now want 
to analyse the system (|4|,[|) in second order. Let us first neglect the time dependence of the 
coefficients and make a plane wave ansatz 

V oc M oc e^ k • . 

Inserting this ansatz into (HJB]), we obtain the following dispersion relation in lowest order Ht 
and ujtx, i-e., for short waves; 

k kH T (R? + {A/h){R + l)) 
w^ + vy , "o= V3{1+R) , 7 = ^ , (6) 

with R = Sps/^Pr @> El- I n the matter dominated regime, R » 1, we have 7 « « ^J 21 
where Z denotes the wavelength of the perturbation. 
We then approximate the total damping by 



with / = I j(t)dt , 



where t en d is defined by tT(t en d) = l/2ir or £end = tdec (whatever condition is satisfied earlier); 
and ti n is the time when the perturbation enters the horizon ti n ~ 1/2 or the time of reionization, 
tin = U (whatever happens later). 

If we use this estimates for damping of fluctuations induced by a spherically symmetric 
collapsing texture [||, [5], g|, we obtain 



f 1+2^ 3/2 



V 1 + J 



for z c > z de 



else, 



(7) 



where z c denotes the redshift of texture collapse, t(z c ) = t c I and we have assumed z% > z c . 
This naive estimate, which leads to a factor 15 damping for z c = 90zrf ec (t c = 0.1^ ec ); 
a factor 5 damping for z c = 10% ec (i c = 0.3trf ec ) and a factor 1.6 of damping for 
z c = 2z dec (t c = (l/V2)t dec ). 

In Fig. 1-4 we show the results from numerical integration of the system (|],[|), with a 
gravitational field induced from a spherically symmetric, collapsing texture [10|. Comparisons 
with the naive predictions (0) show that our approximation is actually quite reasonable! 

This result tells us also that we must require at least z% > 2z^ ec to obtain damping by, say, 
a factor of 2. 



2.2 Spectral distortion 

In addition to damping the amplitude of fluctuations, the reionized plasma also distorts CMB 
spectrum: The ionizing radiation does not only ionize matter but also heats up the plasma 
to temperatures around typically T e ~ 0.3 to a few eV. Nonrelativistic Compton scattering 
up scatters the low energy microwave photons. If the Plasma has a thermal (Boltzmann) 
distribution and the Photons are Planck distributed (as in our situation), the induced change 
in the spectrum can be described by a single parameter, the Compton-y parameter 12, 13]: 

y = — n e (T e - T CMB )adt « 0.4 x lO" 5 — r —*- (T e /2eV) , (8) 
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where T e is a weighted "mean electron Temperature" : 

T e = ( Zi + I)- 3 / 2 H T e {z + l) l ' 2 dz . 
Jo 

(This approximation is excellent for T e >> Tcmb and Zj >> 1, i.e. the situation we are 
interested in.) 

The observational limit set by the FIRAS experiment on COBE (see J. Mather, this Pro- 
ceedings) is 

y < 2.5~ 5 . (9) 

In the next section we shall see, that ionization can be maintained by collisions only if T e > 1.5eV 
and, on the other hand, it seems to require quite some fine tuning to maintain the plasma ionized 
by photoionization without heating it up to a temperature of about (1 - 2)eV. If we therefore 
assume a mean electron temperature of 

T e > 1.5eV , 

observation @ already limits the ionizing redshift to 

Improving the observational limit of the y-parameter by about a factor of 5 would there- 
fore rule out reionization which happens early enough to lead to significant damping of CMB 
fluctuations {z% > 2zd e c, say)! 

3 Requirements for an ionizing radiation field 

We now assume that at some high redshift z >> z^ ec there exists an ionizing radiation field 
which originates from the first nonlinear density perturbations, i.e. the first 'macroscopic' 
baryonic objects in the universe formed by some unknown mechanism, and which is capable of 
reionizing the universe. We want to study the properties of this radiation. 
Its energy density pi and its spectrum /j are related by 

1 f°° 

Pi = ^ u; 3 fi(uj)dLO (% = c = k B = 1) ■ (10) 

7T Jo 

We set 

Pi = qpb 

and assume that q is slowly varying over some time period, where ps is the baryon energy 
density. For this radiation to be able to ionize the universe, we must require 

ni(uj > A) > n B (A = lRy = a 2 m e /2 = 13.6eV), where (11) 
1 f°° 

Hi = I u) 2 fi{u))du , 
jo 

a denotes the fine structure constant and m e is the electron mass. To be specific, let us assume 
that fi is a Planck spectrum with chemical potential p >> 1, 
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Table 1: Limits on diffuse background radiation from X-rays down to the near infrared are 
listed with the corresponding references. The precipitous rise of the limits below lOeV is due 
to contributions from airglow and zodiacal light 



For this spectrum one finds 

6 6T 4 e _/i 

Pi ~ — ^T, e P and thus q = — | . 

ir z n Pb 

The requirement ([ll]) then yields 

g>0.8x Hr 7 e (A/Tl) [(A/Ti) 3 - 2(A/Ti) 2 + 2(A/T i )]~ 1 



(12) 



This condition is shown by the heavy line in Fig. 5. It is compared with observational limits on 
background radiations at the wavelengths corresponding to the temperatures given in Table 1. 
The lower light line, which corresponds to those limits today, 2 = 0, shows that if there exists 
a unresolved ionizing radiation background today its energy has to be close to about 2eV and 
its intensity is only little (less than a factor of 3) below the observational limit. In the upper 
light line we translate the limits of Table 1 to z = 10 (multiply by 2 + 1). The dashed region 
is the allowed parameter range if the universe is to be photoionized until z = 10, i.e., until 
recombination is slower than expansion and the universe remains ionized without a ionizing 
radiation. 



3.1 Processes in the plasma 

We now calculate the ionization and recombination rates and the heating and cooling functions 
for the cosmic medium, neglecting the Helium contribution: 
Photoionization: The Karzas-Latter Photoionization cross section is given by 

647T 

^ = 2 I U ^/dfgbf , (13) 



m%a?>\pi> 



where g^f is a Gaunt factor which we set equal to 1 in the sequel |TJ] . The photoionization rate 
is correspondingly 

i 8irm e a 5 f 00 _ x dw 



G 



where x = n e /nB is the ionization fraction. Inserting now the thermal spectrum defined before 
we obtain 

t- 1 » 6s" 1 x g(l - x )(A/T l )*E 1 (A/T l )(^)(^f . (14) 
Here E\ denotes the usual integral exponential function 

poo „— x 



roc q—x 

Ei(x) = / dx . 

Jx X 



The photons do not only ionize the plasma, but with their remaining energy they also heat 
it up. This photoionization heating, F p i is given by 

87rm e a 5 ,„ N f°° to - A „ dcu 



« 78(eV/s) x q(l - *)(^)(^)3 (A /T 4 )V A / T * - (A/T^A/T,)] . (15) 

Recombination: We obtain the the recombination cross section from a detailed balance 
argument and the photoionization cross section p9| . This leads to the recombination rate 

_ i ^ A l/2 2 13/ 27r l/2 a / 3Te /2A + l \ 2 

tr " ^33/2 I (r e /A)V2 ) x nB 

« lO- 13 ^ 1 x (A/T e ) 1 /V(^)(i^) 3 . (16) 
v 7 ' K 0.025 A 200 7 v ; 

Due to the loss of one free particle, recombination leads to a loss in kinetic energy, a 
contribution to the cooling function A: 

A r = (SfflTjI 1 = 2 x 10- 12 — (r e /A)VV(^-)(— )3 . (17) 
V / ! r g \ e, ) 0.025 200 ' v ' 

Collision: Collisions of energetic electrons can lead to ionization and/or excitation of 
hydrogen. For high enough electron temperatures this is a very efficient cooling mechanism 
even is only a small fraction of neutral hydrogen is present. The collisional ionization rate is 



about (T|, H| 



t" 1 = 6 x 10-V 1 x (T e /A)VV^(l - *)(^)(^) 3 • (18) 
The collisional ionization cooling rate is 

A c = 8x 10- 7 — x (Te/A^e-A/Te^ _ x)( M!)(i±f )3 . (19) 
s y 1 1 K A 0.025 A 200 ; y ' 

Correspondingly, for the excitation cooling one obtains 

Excitation cooling is thus always nearly 10 3 times faster than ionization cooling. Hence, the 
latter can be neglected. 
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As long as the electron temperature is substantially below leV, excitations are not common 
enough to cool the plasma efficiently. In this situation Compton cooling of the CMB photons 
is the most efficient cooling mechanism. 

Compton cooling: From the Kompaneets equation one obtains 



nB^CMB = ^ m — ( U( f T ) / U 4 fcMB(fcMB + l)duJ 
1 7T z m e Jo 



-10 eV „J T e~ T CBM \ ( 1 + Z 



K ur^-s-^j^j . (21) 

The corresponding heating process by compton scattering of the ionizing radiation can be 
neglected (since there are so much less ionizing photons than CMB photons). 

Of course the plasma is also adiabatically cooled due to expansion and it is important to 
compare this with the cooling rates above. Furthermore, if the process of expansion is much 
faster than one of the processes calculated above, we can neglect the corresponding process in 
the expanding universe. 

Expansion: In a matter dominated Friedmann universe with Q t ot = 1, the expansion rate 
is given by 



H ss 4.6 x KT 15 ^ 1 x (h/0.5) ( ^— -f ) . (22) 



exp - - ■• v-, V 20Q 



l + Z 



3/2 



Adiabatic cooling due to expansion amounts to a cooling rate of 

fdT e \ n ..„„_i a eV , AW1 , ,fl + z\ 3/2 

ad 

An additional possible cooling mechanism would be Bremsstrahlungs cooling but this turns 
out to be very small within the range of electron temperatures we are interested in. 



A exp = (3/2)(l + x) (^) « 2 x 10- 13 ^ x (T e /A)(l + x) {^)" 2 (h/0.5) . (23) 



3.2 The differential equations 

The degree of ionization, x, the electron temperature, T e , and the form of the ionizing spectrum, 
fi are in principle determined by the following system of differential equations: 

^ = -t^+t-'+t- 1 (24) 

<El = 2^T e + - 2 - (r-A)--^-^ (25) 
dt a e 3(1 + x) V ' 1 + x dt v ' 

dfi_adfi_ _ _djj L dfi _ dfc 

dt adu; ~ dt lpt + dt lrec dt \ Rom P aneets + source ■ 

The first term on the right hand side of the second equation is adiabatic cooling due to 
expansion. The second term denotes the other heating and cooling mechanisms: 

r = T pi , A = A r + A e + Acmb ■ 

The last term is due to the increase of independent particles by ionization which lowers the 
kinetic energy per particle and therefore the temperature. The third equation sketches the 
changes in the ionizing spectrum which are relevant at energies u > A. The second term on 
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the left hand side accounts for the redshift of photon energies due to expansion. The first 
term describes the loss of photons due to ionization. Then there comes the gain of photons by 
recombination, the energy losses due to Compton scattering off the lower temperature electrons 
and finally the source term which represents the emission of ionizing radiation from some 
unknown primordial objects. This term has to be guessed. 

We have not yet managed to solve the system of differential equations (24,25,26) in some 
generality. We have just looked at the stable situation, ^ = = ^ = and solved 
the resulting algebraic equations for x and T e for a given ionization spectrum /j which we 
parametrized like in Section 2 as a Planck spectrum with chemical potential, fi is then fully 
determined by its temperature Tj and its chemical potential /x or the parameter q which was 
introduced in Section 2. The resulting electron temperature and the degree of ionization are 
shown as functions of Tj in Figs. 6 and 7. Of course this is not a realistic situation since 
photoionization is so fast that it will quickly lead to a depletion of the spectrum above 13.6eV 
which is probably not readily refilled by a realistic source term. This leads to the very unphysical 
behavior that even for very small Tj, i.e. very few ionizing photons, x ~ 1. 



4 Results, Conclusions 

Our preliminary results show that early reionization is still possible even though the parameter 
space is squeezed by the limit on the Compton y parameter and the absence of any near infrared 
background. The first fact tells us that the electron temperature cannot have been higher than 
2eV for a substantial duration at high redshift, z > Zd ec . The second constraint reveals that, on 
the other hand, there was no ionizing uv-radiation present at z ~ 10 which might have kept the 
intergalactic medium ionized. If is was ionized at this redshift it thus was collisionally ionized, 
i.e., the electron temperature was above 1.5eV. 

To weaken these constraints, one might still imagine the cosmic plasma to be photoionized 
up to a redshift z < Zd ec , say z 50 and collisionally ionized later. But it is not clear if 
even this possibility remains, i.e. if it is possible to constantly photoionize the plasma without 
heating it above 1.5eV, say. To decide on this last possibility we have to fully solve the system 
of differential equations presented in Section 3.2 for some reasonable source models. In addition 
we have to take into account the clumping which must be present at a time when the first formed 
objects are supposed to emit ionizing radiation: Whenever a factor n 2 B enters our equations 
(e.g. recombination and collisions) we have to replace the background Friedmann value of n B 
by 

n (clump) = n | (1 + J 6{x) 2 d 3 x) = n | (1 + J p {k)d 3 k)j 

where 5(x) denotes the density fluctuations in the baryon distribution and P(k) is the power 
spectrum, P(k) = (S(k)) 2 . 

Furthermore, if the primordial plasma is contaminated by metals from the first objects, these 
might contribute to the cooling rate substantially and thus influence the electron temperature. 
It is an important but difficult task to estimate this effect. 
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Figure Captions 



Fig. 1 The induced microwave anisotropy for a texture collapsing at t c = 7.6 is shown in the 
case of a reionized universe (full line) and a not ionized universe (dashed line) as a function of 
the impact time r. In the units chosen, the decoupling time is i^ec = 74, so that in this case 
Zc/zdec ~ 90. The damping in this case is approximately a factor of 15 and the damped signal 
is also widened to a width of approximately t^ec- 



Fig. 2 Like Fig. 1 but for t c = 20, thus z c /zd ec ~ 13. The damping factor is approximately 
5 in this case. 



Fig. 3 Like Fig. 1 but for t c = 40, thus z c j ' z^ ec ~ 3. The damping factor is approximately 2 
in this case. 



Fig. 4 Like Fig. 1 but for t c = 100 > td ec . No damping occurs in this case. (The two curves 
overlay.) 



Fig. 5 The limiting Pi/ps = q is given in units of 10 -7 which leads to enough (> rig) 
ionizing photons to reionize the universe for a given radiation temperature Tj (plotted in units 
of the ionization energy A = 13.6ey, heavy line. This limit is compared with the limits on 
background radiations at different wavelengths given in Table 1 (crosses connected with light 
lines). The lower limit applies if the ionizing radiation was emitted until today, z = 0. The 
upper limit applies if the emission of ionizing radiation seized at z = 10. The remaining dashed 
region gives the allowed range of parameters in the universe was photoionized until z = 10. 



Fig. 6 The plasma temperature T e is given as a function of the temperature of the ioniz- 
ing radiation for q = 10 and /j a blackbody spectrum with chemical potential for z = 200 
and a static solution of the equations. This figure was produced for his diploma thesis by P. 
Ehrismann[22|. 



Fig. 7 The degree of ionization x is given as a function of the temperature of the ionizing 
radiation for q = 10 and /, a blackbody spectrum with chemical potential for z = 200 and a 
static solution of the equations. The value x = 1 for small Tj is unphysical. It is due to the 
staticity assumption. This figure was produced for his diploma thesis by P. Ehrismann [22]. 
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